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Abstract. We consider the electron transport properties through fully interacting 
nanoscale junctions beyond the linear-response regime. We calculate the current 
flowing through an interacting region connected to two interacting leads, with 
interaction crossing at the left and right contacts, by using a non-equilibrium Green's 
functions (NEGF) technique. Tthe total current at one interface (the left one for 
example) is made of several terms which can be regrouped into two sets. The first 
set corresponds to a very generalised Landauer-like current formula with physical 
quantities defined only in the interacting central region and with renormalised lead 
self- energies. The second set characterises inelastic scattering events occuring in 
the left lead. We show how this term can be negligible or even vanish due to the 
pseudo-equilibrium statistical properties of the lead in the thermodynamic limit. The 
expressions for the different Green's functions needed for practical calculations of the 
current are also provided. We determine the constraints imposed by the physical 
condition of current conservation. The corresponding equation imposed on the different 
self-energy quantities arising from the current conservation is derived. We discuss in 
detail its physical interpretation and its relation with previously derived expressions. 
Finally several important key features are discussed in relation to the implementation 
of our formalism for calculations of quantum transport in realistic systems. 



PACS numbers: 73.40.Gk, 71.38.-k, 73.63.-b, 85.65.+h 
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1. Introduction 

Electronic transport through molecular nanojunctions exhibits many important new 
features in comparison with conduction through macroscopic systems such as bulk or 
thin layers of semi-conducting molecular crystals as used in conventional molecular 
electronics. In particular, local interactions, such as Coulomb interactions between the 
electrons and scattering from localized atomic vibrations, become critically important. 
In crude terms, these interactions are more pronounced in nanoscale systems because 
the electronic probability density is concentrated in a small region of space; normal 
screening mechanisms are thus ineffective. Developing a theory for the non-equilibrium 
electronic quantum transport through such fully interacting nanoscale junctions is a 
challenging task, especially when thinking in terms of applications for future nanoscale 
electronics. 

Having a simple expression for the electronic current or for the conductance of a 
nanoscale object connected to terminals is most useful. This is provided by the Landauer 
formula [I] in the form of an appealing intuitive physical picture, which describes the 
current in terms of local properties of a finite region (transmission coefficients) and the 
statistical distribution functions of the electron reservoirs connected to the central region 
C. However, in its original form the Landauer formula deals only with non-interacting 
electrons. This formalism has been used in conjunction with density-functional theory 
(DFT) calculations for realistic nanoscale systems [2l[3llll[5l[6l[71[8l[9l[T0l[IIl[l2]. It 
has helped tremendously for the qualitative understanding of the transport properties 
of such realistic systems, though only on a semi-quantitative level as the calculated 
conductance is often one or two orders of magnitude wrong. The apparent success 
of such an approach relies on the fact that DFT maps the many-electron interacting 
system onto an effective non-interacting single-particle Kohn-Sham Hamiltonian suited 
for the Landauer scattering formalism for transport. However there are many cases 
when such a mapping becomes questionable: for strongly-interacting electron systems, 
low dimensional interacting electron system, strongly coupled electron-phonon systems, 
to cite only a few. 

Furthermore, the single-particle framework of the Landauer approach cannot be 
directly transferred to the many-body context by expecting that a proper inclusion of 
many-body effects in the single-particle energy levels will suffice. In fact, it has been 
shown that there are many-body corrections to the Landauer formula which cannot be 
formulated in terms of single-particle transmission probabilities [131 E] • 

When a single-particle-like scheme is still valid even in the presence of interaction, 
the Landauer approach can be extended to include inelastic effects by using inelastic 
scattering theory in a generalised Fock-hke space [151 UHl El [HI [HI [20l [211 122] . 

In the context of DFT, it has also been shown that the exchange-correlation part of 
the interaction that leads to the presence of an extra Vxc potential is actually introducing 
corrections to the Landauer-like current. These corrections are crucial and need to be 
taken into account when working with DFT calculations [23| |2^ . 
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The Landauer formula has been built upon by Meir and Wingreen [25] to 
extend this formalism to a central scattering region C containing interactions between 
particles, while the left (L) and right (R) leads are still represented by non-interacting 
electron seas. The current is then expressed in terms of non-equilibrium Green's 
functions (NEGF) and self-energies, and in the most general cases it does not bear 
any formal resemblance to the original Landauer formula for the current [261 EHl l27] . 
Other generalizations of Landauer-like approaches to include interactions and inelastic 
scattering have been developed, see for example Refs. [221 EHl EHl [I3] . 

However, in real systems the interaction is present throughout the entire system, 
even at the nanoscale. This is even more true for the long-range Coulomb interaction 
between bare electrons. Hence it is difficult to consider that in the real system all kinds 
of interaction will stay localized or will be sufficiently well described by localized effective 
interaction in the C region only. Another reason is that in time-dependent transport, 
the external field in the leads may not necessarily be screened instantaneously. Transient 
times can be of the same order as the plasma oscillation (a collective many-body mode of 
charge oscillation) period in leads and play an important role in the transient transport 
properties of the nanojunction [SU]. Taking the interaction into the whole system is 
vital. To achieve this, the so-called partition-free scheme has been developed ^1]. It 
allows in principle the calculation of physical dynamical responses and to include the 
interactions between the leads and between the leads and the central region in a quite 
natural way [321 133] . 

In this paper, we use an alternative approach, and generalize the Meir and Wingreen 
formalism (so-called partitioned scheme) to systems where interaction exists in all the 
L, C, R, as well as at the interfaces between the three regions. Since the choice of 
location of these interfaces is purely arbitrary and since the interactions exist at and on 
both sides of the interface, our approach is equivalent to a partition-free scheme. 

However, while keeping the approach and the NEGF formalism of the original work 
of Meir and Wingreen [25] , we derive the most general expression of the current for the 
fully interacting system. From this, we can recover all previously derived transport 
expressions or corrections when introducing the appropriate level of approximation for 
the interaction. We also derive and study in detail the current conservation condition 
and the constraint that it imposes on the interaction self-energies. 

Our formalism also introduces naturally the generalization of the concept of the 
embedding potential when the interaction crosses at the boundaries. With this new 
concept and with the condition of current conservation, we can explore different levels 
of approximation for treating the interaction in the different parts of the system, as 
well as at the interface. We can then check which approximations are more suitable for 
practical numerical calculations of realistic systems. 

Although a preliminary account of our formalism is already given in Ref. |34j, we 
provide here the full detail of the derivation and a much more detailed discussion of the 
physical meaning of our results as well as of the implementation of the present formalism 
for realistic calculations. 
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Figure 1. Schematic representation of a central scattering region C connected to the 
left L and right R electrodes, with respective quantum-state labels {A}, {n}, {p} for the 
three L,C,R subspaces. The electronic coupling of the region C to the L{R) region 
is given by Vlc/CL {yRC/CR)i and the many-body interaction is represent by T,^^ 
within all regions as well as across the LC and CR interfaces. The LC/RC interfaces 
can be chosen to be located at the contacts between the molecule and the leads, or 
inside the leads (as shown above). 



The paper is organised as follows: We start with the description of the system and 
the notation in Section [2j Then we derive the expression for the current for the fully 
interacting system in Section [3] We provide the full derivation of the expressions for 
the different Green's functions needed to calculate the current in Section HI We derive 
the conditions imposed by the constraint of current conservation in Section |5} We 
finally conclude and discuss open questions as well as different schemes to perform the 
calculation for realistic systems in Section [6j We recall some properties of the Green's 
functions and self-energies in [Appendix A|and B, and provide the proof of an important 



relation for the current conservation is given in Appendix C 



2. System model 

The system consists of two electrodes, labelled L and R for left and right respectively, 
which connect a central region C via a set of coupling matrix elements Vlc,rc- The 
interaction — which we specifically leave undefined (e.g. electron-electron or electron- 
phonon)-is assumed to be well described in terms of the single-particle self-energy S'^^ 
and spreads over the entire system. 

For the calculation of the current, we introduce two interfaces LC and RC defining 
the three regions L,C and R and use different labels to name the electronic states on 
each sides of these interfaces. The labels {X, \'}, {n,m}, {p, p'} are used to represent 
the complete and orthogonal set of states for the L, C and R regions respectively. 

In the following, we will use either the full notation or a compact notation for the 
Green's functions G and the self-energy S. Following Figure ([T]), the matrix elements of 
Mij of a Green's function or of a self-energy are annotated Mc for the matrix elements 
of the central region M^^- We use the notation for Mxy and M/j for Mpp/. The 
matrix elements are also annotated Mlc for Mxm or Mql for M„a' and Mrc for 
or McR Mnpi. For the hopping matrix elements we use V\m or Vn\i for Vlc or Vcl at 
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the left interface and Vpm or Vnp' for Vrc or Vcr at the right interface. 
3. The current formula 

3.1. Non- equilibrium Green's functions on the Keldysh contour Ck 
The Green's function on the Keldysh contour Ck is defined as 

G(l,2) = -i(rc,-vl/(l)vl/t(2)), (1) 

where (1, 2) stands for a composite index for space xi^2 (i-e. states X,n or p in the 
L,C or R regions respectively) and time ri 2 on the time-loop contour Ck- The time 
ordering Tc^ of the product of fermion creation (\Ef''') and annihilation quantum 
fields is performed on the time-loop contour Ck [33, ESI EH EH] • 

The Green's function obeys the equation of motion on the contour Ck [36l [38] : 

[tdr, - hoil)]Gil, 1') = 5{1 - 1') + J d3S(l, 3)G(3, 1'), (2) 

and the adjoint equation reads 

- ho{r)]C{l, 1') = 6(1 - 1') + y" d2G(l, 2)E(2, 1'), (3) 

where ho{l) is the non- interacting Hamiltonian. 

3.2. The current 

Using the continuity equation Vj(l) + dtn{l) = 0, one can write down the current 
through the interface between the L and C regions. It should be noted that considering 
an interface between the L (R) and C regions is a merely a virtual partitioning for 
mathematical convenience — the interactions are present throughout the entire system. 
It is also useful to consider such LC/RC interfaces in order to connect our results 
to previously derived expressions within the partitioning scheme. The location of 
these interfaces is also arbitrary, but for convenience in future numerical computation, 
different choices are possible: for example at the contacts between the leads and the 
ends of the molecule, or at the contacts between the leads and the so-called extended 
molecule, which already contains part of the leads. 

After integration of the continuity equation over the (half) space of the L region, 
one obtains the following expression for the current flowing through the left interface 
(from L to C): 

/L(t) = -e^(iV,(t)), (4) 

where {NL{t)) is the total number of electrons in the L region. It is obtained from the 
lesser Green's function as 

(^^W) = E-i^AA(t,t)- (5) 

A 
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To calculate the time derivative of Gf^{t,t), we first go back to the full two-time 
dependence of G"^, and then take the equal-time limit after performing the derivatives: 

|(iV,(*)) = i: {-'^GUtut) - ,±GUtM^ (6) 

Using the equations of motion, Eqs. ^ and (|3]), for the Green's functions on Ck, 
one obtains the current II as 

h{t) = ^Tr, [(SG)<(t,t) - (GS)<(t,t)] (7) 
where we have re-introduced the h to have the correct units of current and conductance. 



Using the rules of analytical continuation (see Appendix B), we get 

Tr,[. . .] = Tr, [(S^^'<G"^) + {{Vlc + S'^^'O 

- {G< {Vlc + S^^''')) (t,t) - (G"-S^i^'<) it,t)] . ^ ^ 

There are no lesser and greater components for Vlc since its time dependence is local, 
i.e.VLcit,t') = VLcimt-t'). 

In the steady state, all double-time quantities X{t,t') depend only on the time 
difference X{t — t'). One obtains the following expression for the current II after Fourier 
transform, 

h = ll ^Tr, [(SMG(u;))< - (GMSM)<] (9) 

or equivalently 



/l = ^ 1 ^Tr, [VlcG<{u) - G<{u)Vlc + S^^'<(a;)G"^(a;) + S^^'^(a;)G 



(10) 



Now that we have sorted out the Keldysh components, we have to sort out the index 
(matrix elements) parts. First we concentrate on TrA [(S^^G)*-] and TrA [(GS^^)"^] . 
We thus have (we use the symbol ^ for summation to have a better graphical distinction 
between the sums and the self-energies S): 



TrA [(S^^G)<] =E(SrG„A)< +J2i^f^G,„)< +J2i^f^G,,)< , ^^^^ 

A,n A, A' A,p 

and similarly for TrA [(GS^^)<]: 

TrA [(GS^^)<] =J2iG.n^^^)< +E(GAA'S^f)< +E(GapSJ1^)<. (12) 

A,n A, A' A,p 

In this present version of the theory, we assume that the matrix elements and S^J^^ 
vanish, as there is no direct interaction between the left and right electrode. Because of 
the geometry and the heterogeneity of the nano junctions and because of the different 
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dimensionality of the leads and the central region, we assume that there is an effective 
screening of the Coulomb interaction so that the electrons of the L and R leads do not 
interact directly via any or T,^^ matrix elements. The distance |xa — Xp| between 
two points in the L and R leads respectively is large enough so that the spatial decay 
of S^J^^ = E^^(|x;v — Xpl) make the contribution of zero or negligible. In other 
words, the presence of the L and R electrodes affect directly the central region C via 
and ^^c- However the L electrode do not affect directly the R electrode (and 
vice- versa), but only indirectly via exchange and correlation effects involving the states 
of the central region C. This assumption seems to be valid when the size of the central 
region is of the order of several atoms (i.e. a molecule). If the central region were to 
be a single atomic impurity coupled to two continuum of delocalised electron state, this 
assumption would not be valid. 



We now turn to the evaluation of the sums in Eqs.(ll) and (12). First we 
concentrate on the ^ sums: 



^MB\< 
'A J 



l^AA'<^A'Aj -f ^(CtAA'^A'; 
A,A' A,A' 

Ev^MB,<x-Ya I sr^MB,r^< ^< sr^MB,a v^MB,< 
^AA' '-'A'A "t" ^AA' "-^A'A "-^AA'^A'A "-^AA'^A'A 

A,A' 

Ev^MB,</^a r^r \ , /sr^MB,r y^MB,a\^< 

^AA' I'-^A'A - (-^A'aJ + l^AA' ~ ^A'A J*-^A'A 

A,A' 



A,A' 



^AA' *-^A'A - ^AA' *-^A'A " -"-r^ 



^MB,>^< _ Y^MB,<Q> 
L L L L 



A,A' 



In the second line, we have used the rules of analytical continuation. In the third, 
we have used the equivalent of cyclic permutation in the calculation of a trace, i.e. 
swapping the index A and A' in the last two terms. This is possible here since the sums 
and all matrix elements are defined in the single subspace of the L electrode. The final 
result looks like the collision terms usually obtained in the derivation of the generalized 
Boltzmann equation from quantum kinetic theory. They correspond to the particle 
production (scattering-in) and absorption or hole production (scattering-out) related 
to inelastic processes (i.e. non-diagonal elements of the self-energy on the time-loop 
contour E^) occurring in the left electrode. It has also been shown that the integration 



of such term vanishes as a result of the gauge degree of freedom [39]. In Section 3.4 
we use another route to show how and why these terms can vanish by using generalised 
non-equilibrium distribution functions. 



Now let's concentrate on the sums in Eq.(ll) and Eq.(12). Once more using 
the rules of analytical continuation, we can see that we need the knowledge of the 
following Green's functions matrix elements: G^^, G^^, C^y^ and For this we use 

the Dyson-like equation defined for the non-diagonal elements: G^;^ = {n\{GT^g)^\\) , 
We shall not go into the detail of the full calculations for all four Green's function 
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matrix elements; instead we concentrate on one matrix element {n\{G'Eg)'^\X) to show 
the mechanism of the derivation: 



n|(GS^)<|A) = (n|G"E</ + G<SV + G"'S"^<|A) 

= 5Z ^nAi^AiA2fi'A2A + ^^np^'pMdXiX + ^nm^mA2fi'A2A 
Ai,A2, p,m 

+ ^nAi^AiA2fi'A2A + Gnp^pX29X2X + ^nrrX'mX29X2X 
+ ^nAi^AiA25'A2A + ^np^pA2fi'A2A + ^nm ^mA2 fi'A^^Ai 



(14) 



with S'^r = Vmx + ^"^^ ^mA ~ ^oia'^- ^ne has to keep in mind that we use a 



''mX 

model such that S^^^ 



0. 



As we show in detail in Section 4.1, the interaction defined within the subspace 
of the L lead YfJ\ can be factorized out and included in the renormalized Green's 

Ai A2 

functions gx{^x^ of the L lead. Hence the matrix elements G^^ — (^l(G'S5f)^|A) can be 
recast as G^^^ = {n\{Gc ^cl with 

^nX = V^^nm ^mA' ^A'A + ^nm ^mX'^ 9x'X + ^nm ^mA2 ^A'A' (^5) 
m,X' 



Similarly, we find that 



^Xn 



{X\{h ^Lc Gc)''\n), 
Gl^={X\{h ^Lc GcYln). 



(16) 



Finally by using the rules of analytical continuation for the above matrix elements, 
we find that the current fiowing through the left L interface Eq.(lO) can be recast as 

du. 



-Tr 
27r " 

+TrA 



G'cTc'^ + G^(tg''^"f 



^MB,>^< 
L/ L 



.MB,<^> 



(T 



C) 



(17) 



where 



T 



c 



LCi 



~ ^CL {9l 9l) "^LC + ^CL 9l ^lc 



:i8) 



and 



~9l^ 



= ^'cL i~9l - 9l) ^tc - ~9t ^Ic 



(19) 
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Eqs. (17), (18) and (19) are the main results of this work for the current formula. The 



current flowing at the left LC interface is given by two traces: the first trace is a 
generalisation of the Meir and Wingreen expression of the current to the cases where 
there are both interaction within the left electrode and crossing at the LC contact. The 
second trace is a term related to inelastic transport effects involving summation over 
the left electrode states/sites only. 



An expression similar to Eq.(17) can be obtained for the current Iji flowing at the 



right RC interface, by swapping the index L ^ R and changing the sign to keep the 
same convention for positive current flowing from the left to right direction. We find 



— Tr 
+Tr, 



C^C +^C\^C 



Gc 



J- c 



c) 



.MB,>^< 



'-'R 



'R 



(20) 



We now comment on the physical meaning and implication of the new Tg> quantities 
and the different traces entering the definition of the current. 



3.3. Relationships between the Tq quantities 



From the definition of Tq\ Eq. (18), we can also define the quantity T^'^ using the 



L,9 



greater components and g^^ such as 



= {gl - ~gl) + J^hr ^Ic- (21) 
It is now easy to show that T'^^^ and T^^^ are related to each other by 



C C 



(T^)t - T^. 



(22) 



This is a very interesting relationship in the sense that the quantities obey a 



relation of the type: (greater) — (lesser) = (retarded) — (advanced) as for conventional 
self-energies or Green's functions, though T^'^ and T^'^ are not proper lesser and 
greater quantities. It should also be noted that the different T^(ci;) play a similar 
role as the the lead self-energies (or embedding potentials |10l El Hlj), defined as 
S^'^(a;) = VcL gli^) ^Lc when the interactions are not crossing at the contacts. 
However, they are not simply related to the straightforward generalization of these 
embedding potentials. The latter have the following form (see also the expression for 
the Green's functions given in Section 111) 



(Scl(w) 5'l(w) Slc(c^))' 



(23) 



with S 



lc/clK^) 



Vlc/cl + S^(?/c'Li^) 



and X 



<,r,a 



continuation for Y^^'^^ and Y^'"^ do not give the same expression as for T^'' or 



The rules of analytical 
For 

example Y^'"' = ^^l 9l ^lc ^ and ^ tj''. 

This leads us to an important proof of this work: Because of (a) the very existence of 
the interaction crossing at the contact, of (6) the fact that ^1a/ac ^la/ac (as opposed 
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to the non-interacting case where V^^/cfC ~ ^La/ac ~ ^La/ac, and of (c) the rules of 
analytical continuation for products of three quantities, the usual cyclic permutation 
used in the calculation of the trace Tt\[{'EG)'^ — (GS)^] cannot be used to transform 
the initial trace over {A} onto a trace over {n} only. Therefore the current II at the LC 
contact cannot be expressed simply in terms of the generalized embedding potentials 
Yq'^, hence the introduction of the quantities. The current is not obtained from a 
straightforward generalisation of the Meir and Wingreen formula using the embedding 
potentials Yq'^: 

^'^ij £Tr„[r^"'<G> - Yl^'>G<] + Tr,[Sr>G< - Sr<G>]. (24) 
3.4- Non- equilibrium distribution functions in the leads 

Now we consider the terms TrAp^^'^G^; — S^^'^G^], and show in which conditions 
this trace vanishes. For this we introduce the non-equilibrium distribution functions 
fiiuj) and /^^'^(cj) defined from the generalised Kadanoff-Baym ansatz [33] as follows: 

G<{u) = f^{u)Gl{u) - Gl{u)f<{u) 

and similarly for and /^^'^, obtained from G^ and E^^'^. These distribution 
functions follow the conditions fi=fi—^L and /™*'^ = fi^''^ ^ 1l (where the identity 
matrix in the L region is 1l = 5aa' so that the usual relationships between the different 
Green's functions (and self-energies ) X"^ — X" = — X^ still hold. 

For convenience, we consider for the moment that the non-equilibrium distribution 
functions are diagonal in the corresponding subspace, i.e. f^^ = fx^xy- However there 
is no formal difficulty to deal with a full, non-diagonal, dependence of these density- 
matrix-like distribution functions. 



Introducing the definition of the non-equilibrium distribution functions in Eq.(13) 
we end up, after lengthy (but rather trivial) calculations, with 



Tr, 



^MB,>^< _ ^MB,<^> 
L L L L 



(27r)^TrA [(/< - /r'<)Af (a;)Af (c.)J (26) 
where the respective spectral functions are obtained from 

2mAf{u:)=Xl{uj)-Xl{u), (27) 

with X{uj) = G{uj) or S(u;). 



Eq.(26) shows that the collision term is a measure of the deviation between the 
two distribution functions /^(w) and /™*''^(a;). At equilibrium, because all distribution 
functions are equal to the Fermi distribution = f^^'^ = Z*"^, this collision term 
vanishes. At non-equilibrium this is not generally the case. 

One can now imagine the following case: the indices A represent a spatial location 
or a localised electronic state on a lattice point in the left electrode. For A located well 
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inside the electrode, the system is in its local equilibrium (the Fermi distribution with 
a Fermi level shifted by the left bias) and both distribution functions and f^^'^ are 
equal to the left Fermi distribution, hence these indexes do not contribute to the trace. 
Only the lattice points (or states) that are not in local quasi-equilibrium, i.e. those close 
enough to the central region to experience the potential drop and the interaction effects 
at the contact and beyond will contribute to the trace. From a computational point of 
view, this is good news in the sense that one is not be obliged to perform the summation 
in TrA[...] over all the infinite A indexes of the semi- infinite left electrode. Only the 
states/sites for which — f^^'^ 7^ will contribute to the trace. So if we choose the 
location of the LC interface far enough inside the left electrode of the real system, then 
TrA[...] = 0. In this sense, we have just provided a formal proof of the concept of the 
so-called extended molecule [HI [TUl IHl El SSI SS] The extended molecule represents the 
central region C and consists of the molecule itself but also a part of the left and right 
electrodes to which the molecule is connected. The concept of the extended molecule 
has been introduced empirically in realistic calculations of molecular junctions in order 
to avoid any problems with the asymptotic behaviour of the electrostatic potential in 
the bias of finite (not small) applied bias. 

The value of TrA[...] can also be understood as a measure of the "efficiency" of the 
location of the LC interface in the L electrode. The larger (smaller) the value is, the 
farther (closer) from local equilibrium the LC interface is. This measure can help in 
finding a good compromise between having a sufficiently large extended molecule that 
is nonetheless small enough for tractable numerical calculations. 



3. 5. The Meir and Wingreen current formula 

Using different approximations (for example single-particle approaches, mean-field 
theories, interactions localised in the central region only) for the different interacting 



self-energies, we can recover from Eq. (17) all the previously derived current expressions 



for non-equilibrium nanojunctions. We have analysed all these connections in detail in 
Ref. [33] and we will not repeat the analysis here. 



However, as it will be useful below, we now briefiy recall that Eq. (17) bears some 



resemblance to the current expression derived by Meir and Wingreen [25] in the case of 
interaction present in the central region only: 

ie /" dcj 



7-MW 



h / 27r 



Tr„ [/i(G^-G^)r^ + G<r^] 



28 



The connection becomes more apparent when we use the definitions I/lF^ 



VcL 9t Vlc = Sc'^ = -(Sc'<)^ and iF^ = Vc^gl - qDVlc 



jMW becomes 



rMW 



2^ 



-Tr. 



G 



E^;'< + C^r 



c^c 



.L,<\t 



+ G<(S 



and Tp' play the role of the L lead self-energies S 



L,a 
C 



Hence 



(29) 



One can see by comparing Eq. (|17[) and Eq. (29) that the quantities T^, (T^)^' 



, S^''' and E^'^ respectively in the 
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cases where the interactions cross at the L interface. However, as we have discussed 
above, the quantities are the forward generahsation of the embedding potentials 
Y^'^. Note that, in the model of Meir and Wingreen, the leads are non-interacting. 



hence the second trace TrA [...] in Eq. (17) vanishes. 



4. The different Green's functions needed 

For the evaluation of the currents II and I a, we need to know the following Green's 
functions : G"^^'^ and G"j[^^ . For this, we calculate the matrix elements: = 
{n\ig + gJ:Gr\m), Gl = {X\{g + gEGnX'), = {p\{g + gEGrip'). The results 
are given in the following subsections. 

4.I. Renormalisation of the Green's functions in the electrodes 

First we show that the terms in Yf/^^ can be factorized out and included within the 
renormalization of the left lead Green's functions g°[^^'^ ■ Hence the matrix elements 
^nx — ('^K^^fl')^!'^) can be recast as G^^^ = {n\{Gc ^cl with gi the 

renormalised L lead Green's function whose definition is given below. 
We start from 

G<, = {n\{GJ:g)<\X) 

= ^riXi ^AiA2 9X2X + ^nm ^mA2 9X2X + ^nXi ^AiA2 9X2X + ^nm ^mA2 fi'A2A 

Ai,A2,»n 

+ ^nAi ^AiA2 fi'A^.A + ^nm ^mA2 S'a^.A; 



hence (we now use Einstein notation for the sums) 

^nX ^nXi_ ^AiA2 yA2A — ^nXi_ ^AiA2 yA2A "r ^nm ^mX2 yX2X "r ^nm ^mX2 yX2X 

+^nAi ^AiA2 fi'AjA + ^nm ^mA2 ^'AjA? 



(30) 



(31) 



and 



so 



^nAi(l 9"')xiX — C^Ai ^AiA2 S'AaA + ^nm ^mA2 fi'AaA 

+^nm ^mA2 5'A2A + ^nXi ^AiA2 5'A2A + ^nni ^mX2 9X2X^ 

^nX — ^nXi ^XiX2 9X2X ^nm ^mX2 9X2X '^nm ^mX2 9X2X 
+^nAi ^AiA2 fi'A2Ai(l ~ fi'")AiA + ^nm ^mAa fi'A2Ai(l ~ 9°") xl) 



(32) 



(33) 



/AiA' 

where we define the renormalised Green's functions gl for the left L electrode as 

^,^,(1-S^^'V)a'a,=^aa,- (34) 



To solve for Eq.(33) we also need to know the matrix element GJ^^i 



Gl^, = {n\{GT.gr\\^) = (n|G"-S'/|Ai) = G;,^ S^^,^ g^^^^ + G^^ ll^,^.^ g^^,^, (35) 
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hence 



^nXsi^ 5'l)a3Ai — CJ^m ^mAs dx^X^: 



SO 



^nAi ^nm mX2 SX2X1I 

with a similar definition for the renormahsed Green's function ^£ as given 



Eq.(34). With our compact notation, we have 



~r/a 

9l 



r/a . r/a Y^fAliyr /a ~r/a 

9l ■ 



r a , r a ^MJ 

9l +9l 
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(36) 

(37) 
for g1 in 

(38) 



Using the resuh of Eq.(37) into Eq.(33) and using the fact that (1 — 
(1 + El gl), we find that 

^nX = ^nm ^mA2 ^Lai ^AiA2 ^A2A + ^nm ^mA2 ^A2A + ^nm ^mA2 ^A2A 



9l) 



+ <^nm ^mA2 ^^Ai ^AiAz fi'A2Ai(l + ^")aiA + G'^^ ^mX2 9x2Xi{^ + fi'")AiA 



r^r spr ~< , ~a , ^< ^pa ~a 

'^nm ^mX2 9X2X '^nm ^mX2 9X2X ^ nm ^mX2 9X2 

{n\{Gc EcL ^l)<|A), 



(39) 



with 



~9t 



9l 9l + {^ + 9l (1 + 9l)- (40) 

In other words: the are the Green's function in the left L electrode renormalised 



by the many-body self-energy Y^^'^ defined in the same subspace of the left L electrode. 
Similar results can be derived for the Green's functions of the right R electrode. 



4-2. The retarded Green's function Gq in the central region 



We find for G] 



c 



{n\G^\m) 



far - sr-^ - Yc 



L+R,r 



G'c = 9h + 9h Gh + gh Yc 

\9hV - Y^ 



c 



(41) 



where y^^^'^ is the sum = y^'^ -|- y^'^ of the generalised lead self-energies 

Y^''' {a = L,R) defined previously as Y^'"' = {T.ca9a^acY = ^hJL^aC with 

^Ca = Vca + ^Ca • 

Similar expression can be derived for the advanced Green's function Gq in the 
central region by swapping r -H- a. 
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4-3. The lesser Green's function Gq in the central region 
We find for = {n\G<\m) that 

G] 
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G^ 



with Y, 



L+R,< 
C ' 



UJ 



c 



c 



ja=L,R 



(42) 



[T,ca{(^)ga{(^)^ac{(^))'^ ■ Using the 



rules of analytical continuation for products given in Appendix B, one gets Y(j 



L+R,< 



4.4- The retarded Green's function G^j^ in the L and R lead 
We find for = (pIG^p') that 



PP' 



9pp' + f 



-yC L,r f^r 

PPi Pip2 P2P' 



(43) 

where Y^^'^ is an embedding potential of the effects of the central region G (connected 
to the left lead) on the right lead. It is defined as 



ml 



(44) 



with E^(^ 



Vrc + ^/jc'*" similarly for T^^r- -^^^ [[9ci^)]~^ ~ '^c'^i^)]^^ ^ 
retarded Green's function of the central region renormalized by the interaction inside 
G and by the lead self-energy / embedding potential Y^'^ of the left lead only, with 



defined in Eq.(41) as (g^) ^ = {g, 



Similarly we can find the expressions for the L lead Green's function : 



Gl = ~gl + ~glY^''^'Gl, 



Y, 



CR,r 



{9l 



c) 



Yr 



R,r 
C 



(45) 



Finally all the expressions given in this section hold for the advanced G1 ^ by 
swapping r -H- a. 

4-5. The lesser Green's function G^ ^ in the leads 

We first concentrate on = {GUg + g)^^- The calculations are rather lengthy, but 
somehow trival, and include the derivation of many intermediate Green's functions and 
many re-factorizations. We find, in agreement with the results of the previous section, 
that 



and 



G-i 



^CR,x 
~R,r/a 

9c 



1 + GlFf ^O^L (1 + lt"''^G2) + GlY, 



CR,a^a - 



rCR,< 



Gf 



Li 



i^LC9c'^CLy, 



{I 



Y, 



R,r/a 



C 



9c' 



;i + ~9r%^'r9^{i + Ya^'iSn + ~9S''yc'^~9S" 



(46) 



(47) 
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The expression for can be obtained from Eq. (46) by swapping the index L to 



R and the self-energy yj^^'^ by y^^'^ given in Eq. ( |44| 
5. Current conservation condition 

One of the most important physical properties that our formalism should obey is the 
current conservation condition. Before deriving the current conservation condition 
for the fully interacting system, we briefly recall the equivalent condition when the 
interactions are localised in the central region only. 

5.1. Current conservation condition for the case of interaction only in the central 
region 

Within the partitioned scheme devised by Meir and Wingreen, i.e. interaction only in 
the central region C, it can be shown that the following trace 

Tr„ [J:<{u)G>{u) - S>(a;)G<(a;)] = (48) 

vanishes for each u, where ^(0;) is the total self-energy of the region C: S = 
+ + Eg. The a-lead's self-energy is Eg = VcagaVc 



aC- 



Eq.(48) is derived from the definition G^'^(a;) = G''\uj)T?'''^(uj)G"'{uj) in the region 
C, and hence {G'')-\G> - G<)(G"^)-i = S> - S< = - = {G"")-^ - {G')-^ 



With Eq.(48), we can derive a condition that must be fulfilled by the interaction 
self-energy [17] in order to satisfy the current conservation condition + = |17] . 
This condition is given by: 

j duj Tr„ [S^^< G> - S^^> G<\ = 0, (49) 

which means that the integrated collision term must vanish. This is a condition 
familiarly obtained from a Boltzmann-like treatment of scattering theory [47J. When 
the interaction self-energy is derived from the so-called ^-derivable approximation 



it automatically satisfies the condition given by Eq.(49). 



Eq.(49) can also be used as a measure of the accuracy of numerical schemes used 
to calculate approximately the interaction self-energy . It can also be used as a general 
constraint equation in the determination a new functional forms for the interaction 
self-energy. 

Now, when the interaction exists throughout the entire system, the derivation 
described above no longer holds and needs to be generalised to the presence of interaction 
within the leads and crossing at the left and right contacts. 

5.2. Current conservation condition for the case of interaction everywhere 

First we consider the general definition of the lesser and greater Green's functions : 

G< = (1 + G'TI')g<{l + G") + G' S< G". (50) 
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The first term represents the initial conditions before the interaction and the coupling 
between the different regions are applied. 

For the central region, we have chosen the initial condition such as {n\g^\m) = 



(see Eq.(42)). We could have chosen another initial condition. Such choices have no 
effects on the steady-state regime when a steady current flow through the central region, 
however the initial conditions play an important role in the transient behaviour of the 
current [Ml ED E21 ES] • 

For the definition of the lesser left- and right-lead Green's functions, however, it 
is not possible to neglect the initial conditions (before full interactions and coupling to 
the region central are taken into account). It would not be physically correct to ignore 
the presence of the left and right Fermi seas since they are the thermodynamical limit 
of the two semi-infinite leads, and act as electron emitter and collector in our model of 
a device. 



By using the standard Dyson equations for we can recast Eq.(50) as 

G< = C {{g')-^g<{g^y^ + S<) G" = G'' S< G^, (51) 

with = + 7^ and 7^ = {g^)^^g^{g"')~^ , and similarly for G^ . Hence 
7< - 7> = (^'*)-i - (/)-! and S< - S> = {G")-^ - {G'')-\ From these properties, it 
can be easily shown that 

Tr,u [S<G> - J:>G<] = 0, (52) 

for each u. The trace runs over all indexes in the system all = {X,n,p} and the 
interaction S are spread over the whole L, C, R regions. This is a generalisation of 



Eq.(|48j). 

Because the trace runs over all the three subspaces, we can apply the usual cyclic 
permutation and recast Eq.([52|) as follows 



- Tr,u [(SG)< - (GS)<] + Tr,u b^G> - ^>G<] = (53) 

or equivalently 

j du Trail [(SG)< - (GS)<] + Tr^u [7>G< - j<G>] = (54) 

Expanding the trace in the first term over each subspace Traii[...] = Tr^f..] + Tr„[...] + 
Trp [...], one can easily identify the definition of the left II and right Ir currents from 
TrA[...] and Trp[...] respectively (see Eq.([9]) above). 

Hence the condition of current conservation II + Ir = leads to 



da; Tr„ [(SG')< - (G'S)<] + Tr^u b>G< - ^<G>] = 0. (55) 

After further manipulation of the trace Tr„[...] using the rules of analytical 
continuation and the relationship between the different Green's functions and self- 
energies, we find that the current conservation implies that 



da; Tr„ 



.^MB ^ Y^+RyG< - (S^^ + F^^+^)<G>J + Tt,,r b>G< - ^<G>] = 0, 

(56) 
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where in the second trace the sum runs only over the left and right subspaces, because 
we have chosen the initial condition for the central region C such that 7>'< = 0. 



What is really interesting with the first trace Trn[...] in Eq.(56) is that it has exactly 



the same form as Eq.(52), but with all quantities (as well as the trace) defined within 
the central region C only (i.e. S = + l^"*"^ in the region C). By looking at the 
definition of G^'^ given by Eq.(42), we can establish that, equivalently to Eq.(52), the 
trace actually vanishes: 



Tr 



(s^^ + y^+^)>G< 



c 



0, 



for each uj. Therefore Eq.(56) reduces to 



(57) 



(5J 



a condition which however is almost systematically verified (see Appendix C). 



Hence it is better to consider Eq.(57) to find the condition imposed by the current 



conservation. Indeed by treating each contribution in Eq.(57) separately and by 



integrating over the energy, we can introduce the definition of the currents Ii and 
Ir such as 



j Tr„ [S^^>G< - S^^<G>] - il{uj) + A 



«l(w) - iRiuj) + ^iniuj) 



0, 



(59) 



where = e/h J ia{ou)dijj, and Ai^ 



[CO 



+ «l(w). 



Hence the condition of current conservation 1^ + Ir = leads to the final important 
result of this section: 



Tr 



+ Tr„ [{L o R}] 
= 0. 



r) 



'^c^)Gc - {Yc^ - ^c^)Gc + (Tc'' + (T^'')"f)G" 



(60) 



Using the properties of the quantities (see Section 3.3), we can rewrite the 
condition imposed by the current conservation as 



/ duj TrQ,=L^c',iJ 



+ Tr„ 



(F^'> - T^c^)G^ - (Fc - ^c^)Gl + (TJ-' + (tJ;'')^)G" 



(61) 



+ Tr„ [{L ^ R}] 
= 0. 



To understand fully the conditions of current conservation, we make the following 
observations: 
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Eq. (61 ) is the generalisation of Eq. (49 ) for the systems where the interaction spreads 



throughout. Like Eq.(49), it contains similar terms involving the left L and right 



R lead as well. But it also contains terms arising from the fact that the interaction 
is crossing at the LC and RC contacts. However, the physical interpretation of 



Eq.(61) still corresponds to the fact that the total integrated colhsion terms must 



vanish. 



(ii) For interaction present only within the C region, one can show that T^''+(T^'')"'' = 

Vlc, and likewise 



as well as Y. 



L,< 



c 



T 



L,l 
C 



and Y, 



C 



^ C 



since S^^f 



for the terms involving the RC interface. Furthermore, in that case, there are no 



interactions in the leads ji 



0, and hence one recovers Eq.(49) as expected. 



(iii) 



Eq.(61) is also consistent with the one of the main point made in Section 3.3 the 
current Jq, (a = L or R) cannot be obtained from a trace over the central region 



defined as Tr„[yj- 



after Eq. (59) as Ai 



c 



[CO 



Y^'>G<] 



If it were the case, then the Aiaioo) defined 



conservation condition would imply that Eq.(61) reduces to 



0. And the current 



(62) 



This equation concerns the collision terms within each of the three regions, but is 
not complete because it does not show any constraint on the interaction crossing 
at the LC and RC contacts. 

6. Discussion and conclusion 



In this paper, we have derived a complete and exact expression of the current crossing at 
the LC (or RC) interface (defining the contact between the L (R) lead and the central 
region C) for general systems with interaction both within each L, C, R region and 



crossing at the left and right interfaces. Our result for the current Eq.(17) and Eq.(20) 



is general and obtained under only one approximation: there is no direct exchange and 
correlation effects between the left and right lead; a condition that is physically sound, 
especially for a large-ish central region where the spatial gap between the two electrodes 
is large enough so that the two electrodes interact only indirectly via the central region. 

Our formalism includes all the cases previously studied with any kind of interactions 
present in the central region only [Ml[55l|56l|571[58llMllMll6ll|63[6a[Ml 
ESI ESI Eni E] ■ It also include other classes of problems such as those where interactions 
also exist within the leads but does not cross at the contacts [72]. Our formalism also 
provides a natural way to extend cases where the excitations exist in the leads and could 
cross at the contacts between the central region and the leads [731 E] • 

We now discuss in more detail different open questions that are of importance for 
applications of our formalism to realistic systems. 
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6.1. Location of the interfaces 

There is one arbitrary choice in our derivation: the location of the LC and RC interfaces 
with respect to the physical realistic system. Such locations are somehow arbitrary in 
our formalism but could be conveniently chosen for practical numerical calculations. 
We have already seen that, in some cases, a local quasi-equilibrium is reached within 
the left and right leads. In such cases we get simplified results for the current, since 
the local non-equilibrium distribution functions /^^, are equal to the local Fermi 
distribution in the left or right lead fL,R = fi%- 

6.2. The extended molecule 

Our formalism provides a formal justification of the concept of the extended molecule 
that has been used so far within the conventional partitioned scheme with interaction 
present only in the C region. Being general, our formalism also provides the 
corresponding "correction" terms needed when the interactions cross at the contacts 
and when the contacts are not in their respective local (quasi) equilibrium. 

Our work also justifies the scheme recently used in Ref. [H] where two closed 
surfaces are used to define two sets of LC and RC interfaces around the molecule. 
Within the inner region, the electron-electron interaction is calculated via the GW 
scheme with the many-body self-energy S(r, r') = G{t,t')W{t,t'). While the screened 
Coulomb interaction W = v + vPW and the polarization P are calculated for the 
extended region in order to ensure a better treatment of non-local screening effects. To 
some extent, this is an empirical way of defining the concept of generalised embedding 
potentials that we obtain in a formal manner in our formalism. 

6.3. Generalized embedding potential 

In a broader context, our formalism introduces in a formal manner the generalisation 
of the concept embedding potential to interacting cases. In the case originally studied 
by Meir and Wingreen, where the interaction is present only in the central region, the 
effects of the leads on the Green's functions defined with the subspace of the region C are 
taken into account via the so-called lead self-energies : Sg,(c<;) = Vcaga{<-^)Vac- These 
non-local self-energies are simply a matrix representation of the so-called embedding 
potential originally defined in real space by Inglesfield [lOl [75], CHI SU Il2] . The latter 
can be seen as defining a surface (or two interfaces) around the central region that the 
interaction does not cross. The non-locality of the embedding potentials arises only 
from the Green's functions defined on this surface. 

In our formalism, when the interaction can cross the LC / RC interfaces, we 
obtain in a systematic way a generalisation of the embedding potentials Yq., defined 
as Y^{uj) = S(7a(c<;)^a(a;)Sac(w)- These generalised embedding potentials contain a 
"double" non-locality, in the sense that the part of E 

ac can have a larger spatial 

extent than the hopping matrix elements Vac- Hence Yq defines somehow not a simple 
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surface but a "buffer" zone that is contained between two surfaces whose separation 
is related to the characteristic (spatial decay) length of the interaction self-energies 



6.4- Calculation of the self-energies 

For practical numerical calculations of the current, one needs to choose the form for 
the interaction self-energies. They can be obtained from conventional many-body 
perturbation theory and Feynman diagrammatics, extended onto the Keldysh contour. 
The interaction self-energy can be obtained from the ^-derivable conserving 

approximation [IHl HHl [3HI [39] , and then they should automatically satisfies the condition 
of current conservation. We have given an example of such interaction self-energies in 
Ref. pi] for the case of electron-phonon coupling inside the region C and crossing at the 
LC interface. One could also devise other functional forms for the self-energies, such as 
functional of the charge and spin densities, and or of the current density itself. In these 
cases, one should devise the functionals such that the condition of current conservation 
is indeed fulfilled. 

Finally, one should note that once the functional forms of the self-energies are 
chosen, one can perform the calculations self-consistently. The self-energies in the three 
regions L,C,R and at the interfaces LC and RC are functionals of the other electron 
(and phonon) Green's functions (or other physical quantities related to them such as 
the charge, spin or current density) defined inside the L, C, R regions as well as at the 
two LC and RC interfaces. Hence the Green's functions and self-energies need to be 
determined self-consistently in all the parts of the system in order to get the current of 
a fully many-body interacting nanojunction at non-equilibrium. 



6. 5. A special case for the self-energies 

When modelling the self-energies such that T^'' + (T^'')"'" 



0, we can express the current 



Eq. (17) as 



Tr 

+TrA 



(G^^-G^)Tg'' + G<(T^-(T 



^MB,>^< 
L L 



^MB,<^> 



(63) 



which bear even more resemblance to the Meir and Wingreen result, Eq. (28), and hence 



could be recast as a generalized Landauer-like expression for the current [H]. In such 
cases, the condition of current conservation becomes 



+ Tr, 



c 



(64) 



+ Tr„ [{L ^ R}] 
= 0. 
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Eq. (64) looks just like the sum of the integrated collision terms in the form of 

J dio Tr^ [J:f'>G< - J:f'<G>] = 0, (65) 

where S™* represent some interaction self-energy and the sum /3 is over not only the 
three subspaces L, C and R but also the two interfaces LC and RC. 

By using the rules of analytical continuation and the relationships between the 
different Green's functions , one can find that, in general, the sum T^j^q + (T'^^^)''' is 
given by 

Tc'' + (Tg'')^ = imcL ^ic - m^L ^tc + mtc - ^cl mic- m 

There are two different cases in which the sum T^'' + (T^'')"'' vanishes: 

(i) when the interaction is only present in central region. We have already discussed 
that case in length in the paper. 

(ii) when interaction is instantaneous, i.e. local in time E'^^(r, r') = v^^{t)6{t — r'). 
Hence there are no lesser /greater components of the self-energy S^^'^^ = 0. This 
occurs, as already discussed, in mean-field based and in density-functional-based 
theories for which a single (quasi-)particle description of the system is available. 

Finally, it would be interesting to study and find cases which go beyond mean-field or 
density-functional-based methods, and for which T^'' + (T^'')''' = 0. If such cases exist, 
it would still be possible to analyse their transport properties in terms of a generalized 
Landauer-like approach pH ]. 
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Appendix A. Relationships between Green's functions 

By definition, complex conjugation of the different Green's functions follows the rules: 

G"'(l,2) = (^"(2,1))* 
G^(1,2) = -(G^(2,1))* 

Similar expressions also hold for the self-energies S. 

Furthermore, there exists relationships between the different components of the 
Green's functions (or self-energies) on the Keldysh time-loop contour Ck- They are 
given by: 

X" = X++ - x+- = X-+ - 
= X++ - X-+ = x+- - x-~ 

X-+ - X+~ = X'' - X", 



Non- equilibrium many-body transpori 



22 



with X''i''2(i2) = G'?i'?2(12) or S^i''2(12), and where {i = 1,2) is the composite index 
for space-time location (xj,tj) and rji is the index of the Keldysh time-loop contour 
Ck branch (+ forward time arrow, — backward time arrow) on which the time ti 
is located. The conventional lesser and greater projections are defined respectively 
as = X^^ and X^ = X~~^, and the usual time-ordered (anti-time-ordered) as 
X* = X++ (X* = x— ). 

Appendix B. Rules for analytical continuation 

For the following products P{i){r, t') on the time loop contour Ck, 



we have the following rules for the different components P^^{t,t') on the real time axis: 

(x = r,a, >, <) 



Appendix C. Proof that / du Tr^,^ [l^C< - 7<G'>] = 

From the definition = {g^)~^g^{uj){g"')~^, the quantities 7^ have only matrix 

elements within the a = L,R lead since they involve only the non-interacting Green's 
functions. Consequently we only have to prove that 



We now introduce the non-equilibrium distribution f^^'^^u) and f^'^{u) defined, 
in the a lead subspace, from the generalised Kadanoff-Baym ansatz [43j as follows: 






(C.l) 



G':(u;) = /°^'(c^KM-^?;M/rM, 



(C.2) 



with X =>,<. The distribution function f^"^ is defined for the non- interacting 
and uncoupled lead a, and is given by the Fermi distribution of the lead a at its 
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own equilibrium, 
manipulation as 



Using these definitions, we can reformulate Eq. (C.l) after further 



W Tr,,« [{fl< - /<)A^i?J , 



(C.3) 



where 



BJuj) 



{Gl{uj) - G^(a;))/2vri, 



(C.4) 



In the second equality of Eq. ( |C.3 ) we have used a diagonal representation for the 
non-equilibrium distributions and f^. 



There are several different cases for which Eq. (C.l) vanishes: 



(i) Following the same reasoning as in Section 3.4, when the LC interfaces is located 



well inside the lead L, the corresponding states are in their local equilibrium and 



J a 



J a 



0. 



(ii) Following the definition Ba 

Hl — i^))/27ri = ///tt — )■ 0. Therefore Eq. (C.3) vanishes unless A^{u) is singular 
at some energy. This would correspond to the appearance of localised bound 
states in the a lead induced by the interaction (and/or by the non-equilibrium 
condition). Although such an appearance cannot be ruled out in principle, it seems 
to correspond to pathological cases not relevant for the description of the metallic 
leads used in the experiments. 
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